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Abstract 

n , ■ We study existence, uniqueness, and distributional aspects of generalized solutions to the 
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Cauchy problem for first-order symmetric (or Hermitian) hyperbolic systems of partial differ- 
ential equations with Colombeau generalized functions as coefficients and data. The proofs of 
f^ ■ solvability are based on refined energy estimates on lens-shaped regions with spacelike bound- 

"ti ' aries. We obtain several variants and also partial extensions of previous results in [261 1231 [16] 

and provide aspects accompanying related recent work in [281 1101 12| • 
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t^ ■ 1 Introduction 

C^ ' In this paper we establish existence and uniqueness of a generalized solution to the hyperbolic 

Cauchy problem 



O x-^ . 

dtU + ^A^d^^U + BU = F on(0,T)xE", (1) 

U\t=o = G, (2) 

l^% ' where U, F and G are vectors of length m and A^ and B are m x m-matrices whose components 

are generalized functions in the sense of J.F. Colombeau (cf. [3l |4]). All coefficients and data 
may therefore represent functions or distributions of low regularity. Our main focus as well as the 
essential methods are following up along the lines of the seminal papers [111 US] [53] . 
Problems of the type (II])-@ play a prominent role in models of wave propagation in highly 
heterogeneous media with non-smooth variation of physical properties such as density, sound 
speed etc. For more details on motivations from the natural sciences and for further mathematical 
aspects in the context of the theory of generalized functions we may refer to the papers mentioned 
above as well as to the following series of papers on closely related research [5] [33] [T5] [TH] [Ml [13 [H] • 
Second-order wave equations in a similar mathematical context have been discussed in [TT] [35] [T3] . 
Besides existence and uniqueness of generalized solutions to the Cauchy problem we are also inter- 
ested in the relation of the unique Colombeau solution to more classical and weak or distributional 
solution concepts, if the coefficients are of compatible regularity. The analysis of such questions 
and several convergence results are also going back to earlier investigations in [35] [35] [33] [TH] . 
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The outline of our paper is roughly as follows. After a brief reminder of basic notions from 
Colombeau theory of generalized functions in the following subsection, we devote a section to the 
details of the construction of lens-shaped domains and on energy estimates on such domains with 
explicit expressions for the constants. These estimates are then the essential ingredients in proving 
several variants of existence and uniqueness results in Section 3. More precisely, Theorem 3.1 is 
an extension of the main theorem in [23, p. 98] to the case of complex matrices and relaxes the 
required constancy of the coefficients for large spatial distances to boundedness. Theorem 3.4 is 
based on 5^2 -spaces and gives a result for systems of partial differential operators which is similar 
to [m Theorem 3] for the case of scalar pseudo-differential operators. 

The final section investigates regularity as well as compatibility of Colombeau-type solutions with 
classical and distributional solutions in case the coefficient matrices are sufficiently regular. Propo- 
sition 4.1 is an analog of the compatibility proposition in [23^ p. 99] and [16, Corollary 5], whereas 
Proposition 4.2 is a 0°°-va.ria.nt of the regularity result [HI Proposition 6]. In Proposition 4.4 we 
establish convergence of the generalized solution to a weak solution for arbitrary Lipschitz con- 
tinuous coefficients, thereby accompanying the case study with discontinuous coefficients in the 
acoustic transmission problem carried out in [211 Theorem 2.4 and Corollary 2.5]. 
Basic notation and symbols: Let fi C R" be open, 1 < p < oo, k G No, then W'^'P{ft) denotes the 
L^-norm based Sobolev space of order k onil and H''{H) = W''''^{ft). For s G M, the Sobolev space 
if*'(R") is defined by Fourier transform. If F is a Banach space, then C^{[0,T],Y)™' denotes the 
TO-tuples of k times continuously differentiable functions from the interval [0, T] to Y. Similarly, 
L^([0,r], y)™ denotes the m-tuples of square-integrable functions [0,r] -^ Y (in the Bochner- 
Lebesgue sense). If i? is a commutative ring with unit, then Mm{R) denotes the ring of square 
matrices of size m over R with unit given by the identity matrix I^. For any A S iVfm(C), the 
expression ||^||op denotes the operator norm of A as linear map acting on C™. We use ( • , • ) to 
denote the standard scalar product on C™ and || • || for the Euclidean norm. 

1.1 Colombeau algebras of generalized functions 

This section serves to gather some basic notions from Colombeau theory of generalized functions. 
We adopt the topological viewpoint of the construction of generalized functions based on a lo- 
cally convex vector space, developed in [5]. For a comprehensive introduction to the theory of 
Colombeau algebras we refer to [13]. 

Let £■ be a locally convex topological vector space whose topology is given by the family of semi- 
norms {pjljej. The elements of 

Me :- {{ue)e e £;'"'^1 ■.\/i eJ3N e No Pj{u,) ^ 0(e~^) as e ^ 0} 

and 

Me := {{u,), e ^(°'i] -.-ijeJ-iqe No Pj{ue) = 0{e'^) as e ^ 0} 

are called E-moderate and E-negligible, respectively. Defining operations componentwise turns 
Me into a vector subspace of A^^;. We define the generalized functions based on E as the quotient 
Qe '■= M.eI-N'e- If fi' is a differential algebra, then Me is an ideal in M.e and Qe is a differential 
algebra as well, called the Colombeau algebra based on E. 

Let Vt be an open subset of R". By choosing E = C°°{il) with the topology of uniform convergence 
of all derivatives one obtains the so-called special Colombeau algebra Qc°={n) = Q{^)- In the 
current article we will also use the space E = H^{n) = {h e C^{Tl) : d^h e L'^{n) \/a e NJf} 
with the family of semi-norms 

mHHn) = (.Il \\d"h\\l.^n)f' (fc € No), 

\a\<k 

as well as £^ = W°°'°°{n) ^ {h e C°°(fl) : d°'h e L°°{n) Va S N"} with the family of semi-norms 

IMw^^in) = max \\d°'h\\^^,^. {k e No), 



and E = C°°{I x R"), where / is an open interval, equipped with semi- norms 

ll^ll™,K= max WH^^^.^i,) (meNo,A'CCM"). 

To avoid overloaded subscripts we use notations as in \JJ_^ and denote 

Ql-2{^) := GH^(n), Gl-^{^) ■■= Gw°°.'=°(n) and 0(1 x M") := Gf^^^^j^^^y 

Colombeau algebras contain the distributions as a linear subspace. Their elements are equivalence 
classes of nets of smooth functions, G{^) 3 u = [{ue)e]- We say that a Colombeau function 
u is associated with a distribution w e I?'(r2) if some (and hence every) representative (Me)e 
converges to w in T)'{Q). The distribution w represents the macroscopic behavior of u and is 
called the distributional shadow of u. Not every element of a Colombeau algebra is associated 
with a distribution. 

In [27], the subalgebra G°°i^) of regular generalized functions in G{^) was introduced to develop 
an intrinsic regularity theory in t/(n). The subalgebra Ge' of a Colombeau algebra Ge is obtained 
by demanding that the inverse e-power N in the moderateness estimates can be chosen uniformly 
over all derivatives (cf. Definition 25.1, Chapter VII in '27'). For instance, an element u = 
[iue)e] e Gi^) belongs to G°°ii^) if and only if 

The subalgebra G°°{^) plays the same role within G{^) as C°°(J7) does within 'D'{n) and satisfies 
the important compatibility relation 

G°"in)r]V'{n) ^c°°{n). 
2 Standard lenses and basic energy estimates 

A central element of our proof of unique solvability of the Cauchy problem will be L^-estimates 
performed on lens-shaped subsets of the strip [0,T] x M". Similar constructions have been used, 
e.g., in [2], Part I, Section 2.2, in [TS], Section 4.3, and in [7], Section 4.4. Since we are working 
in a generalized functions setting, it is essential to have precise information on all dependencies of 
constants involved in these estimates. For this reason we devote this section to the construction 
of a special variant of lens-shaped domains and to some basic estimates for these types of lenses. 

Definition 2.1. A standard lens C of thickness T > and radii < Ri < R2 is the image set of 
the map ip : [0, 1] x Br^ -^ R"+^ , 

my)^j (erffM,,) for \y\>R, 

where Bj^.^ is the closed ball of radius R2, centered at the origin. We introduce slices of a lens, 
"He '■— i^i^^Bji^), as well as partial lenses Cq = 1Jq<.^<0 "Hr for Q e (0,1]- The latter are 
compact convex subsets of [0,T] x M" with Lipschitz continuous boundary dCs = Hq U T-Lq. 

The standard lens map tp introduced in Definition 12.11 as well as its restrictions ips : Bn^ — )> Jie, 
y H> ip{Q,y) are Lipschitz continuous, but not differentiable at points (6,y) with \y\ = Ri. 
However, since the collection of these points is of Lebesgue measure zero with respect to 'ilJe{BR2) = 
He, they can be ignored when using the lens map to transform integrals over C otT-Lq into integrals 
over the cylinder [0,T] x Br^ or Br^ respectively (cf. [30], Lemma 7.25 and 7.26). Smooth slices 
He are possible by an easy modification of the lens map, but not necessary for our considerations. 

Lemma 2.2. If u £ C{[0,T] x R") and C is a standard lens of thickness T, 

f\u\dVn+l<T sup f \u{Q,-)\dVn, < ^ f \u\ dVn+1 < T f {uiQ , ■){ dVr, . 



Proof. First note that the map ipQ : Bjf^ — > "He is a (global) parametrization of the slice Hq. 
The volume density on Uq is pe(y) = \/det(L'V'e(2/)'^£'^e(y)) = v^det(Dj^f/'(0, vVDyipiQ, y)), 
where Dyip is obtained from the Jacobian Dip by removing the first column, more precisely Dip — 
{Dqi/j Dyi/j). Hence we have 

\detDij{e,y)\ < pe{y)\\De^iQ,y)\\<Pe{y)T. (3) 

With the help of the transformation formula for integrals and using ^ we estimate 

1 

\u\dVn+i= I I \uo^{Q,y)\\detDi;{Q,y)\dydQ 

Br2 

1 

<T f f \uo^Jj{Q,y)\pe{y)dydQ<T sup f \uiQ,-)\dV„, 
J J eG[o,ii J 

Br^ -He 

which proves the first inequality. For the term -js jr \u\ dVn+i we find 

e 

— / \u\dVn+i^— / / \uo%p{e,y)\\de%D'4}{e,y)\dyd£ 

Cb Br^ 

<T J \uoij{Q,y)\pe{y)dy^T J \uie,-)\dVn. 

Br2 He 

D 

For convenience of the reader we also give a full proof of a Gronwall-type estimate in [2] , Appendix 
A, Lemma A. 3, focusing on explicit expressions for all constants appearing in the calculation. 

Lemma 2.3. Let C be a standard lens of thickness T and suppose that u and f are functions of 
class C([0,r] X M") such that for all 6 G (0, 1] 

i / |u(e, ■)\dV,, < J \u{0, ■)\dVn + aJ \u\ dVn+i + J I/I dV„+i (4) 

«e «o Ce C 

with a > 0. Then we have with C := 2Ta, 

\j\u{Q,-)\dVn<e^'^(^j\u{Q,-)\dVn + j\f\dVn+i) VG G [0, 1] . (5) 

We -Ho C 

Proof. We introduce v € C([0, 1]), v{Q) := ^ J^^ \u\dVn+i for 6 G (0,1] and v{0) := 0. In 
addition we put a := J^ \u{0,-)\dVn + J„\f\dVn+i. From Lemma 12.21 we know that < 
^^'(0) < i/Hel"l'^K for e e (0,1]. Therefore w'(e) < a + Cv{e), where C :== 2Ta. By 
integrating the last equation we find v{Q) < aQ + C L v{T)dT and Gronwall's lemma yields 

Cv{0) < Ce'^^ J^ ae^'^'^dr < (e*^® - l)a for all 8 G [0, 1]. Expressing the result in terms of u 
and /, we obtain 

a J \u\ dVn+i < (e^® - 1) ( / l"(0' ■)\dVn + J \f\ dVn+i) Ve G (0, 1]. 

Ce -Ho C 

Using assumption (U), we obtain ([S]). D 



To a first-order operator P{t,x;dt,dx) — dt + X]?=i ^'^ (^: 2;)(9xj + B{t,x) we assign its principal 
symbol a(t, x; r, ^) = tI„i + X^fci ^"' (^' ^)^j- If the matrices A-' are Hermitian, then the principal 
symbol is Hermitian for all directions (t, ^) £ R"+^. At a point {t,x) S ]R"+^ one may then 
define the forward cone T{t,x) as the set of all directions (t, ^) where a(t,x]T,^) is a positive 
definite matrix. A hypersurface is called spacelike (with respect to the principal symbol of P) if 
its normal vector is almost everywhere contained in the forward cone. In the following lemma we 
will construct a lens whose individual slices are spacelike hypersurfaces with common boundary 
dne^{iO,x)\\x\^R2}. 

Lemma 2.4. Let (A-')i<,;<„ be Hermitian matrices such that \\Aj{t,x)\\ < C for \x\ > Ra- 
Then a standard lens C of thickness T > with radii Ri > Ra and R2 > i?i + r(l + 2^/nC) has a 
unit normal vector field vq on Hq (pointing outwards with respect to Cq) satisfying the inequality 

(7;,a(t,a;;z/e(i,^))r/)>ihp V(f,x) G He VO € (0, 1] Vt? e R™. (6) 

Proof. For a lens of thickness T and radii i?i, R2 the normal vector field on T-Lq with 8 £ (0, 1] is 
simply vq = (1,0)-^ for |a;| < Ri and 

Mt.x) = ^ ( ^n^i^' ^ for Ri < \x\ < R2. (7) 

The inequality in ^ is obviously satisfied whenever |x| < i?i. For \x\ > Ri > Ra wc find by 
virtue of (O that 

^^ {R.2~Ri-TC^) > 1 .^^i^encvcr R2> Ri+T(l + 2J^C). U 

A first-order partial differential operator P — ^t+X^iLi ^''^xj +B with smooth coefficient matrices 
is called symmetric hyperbolic if the matrices A' and B arc uniformly bounded together with all 
their derivatives and the principal coefficients A^ are Hermitian. Preparatory for applications to 
Colombeau theory we perform energy estimates for symmetric hyperbolic operators on standard 
lenses. It is important to keep explicit expressions for all constants involved to have precise 
information on their e-dependence in a generalized setting later on. We provide L^-estimates in 
two versions, the second of which can be interpreted as the limiting case for lenses with infinite 
radius. 

Lemma 2.5. Let a symmetric hyperbolic partial differential operator P — dt + '^j^i A^ dxj + B 
be given, where A^ and B are matrices of dimension m. 

i) Let C C [0, T] X R" be a standard lens of thickness T that satisfies inequality ^ and put 
a{L) := 1 + IJdivA — B — -B*||^oo(£\, where divA = X]?=i ^x^AK Here and in the sequel, 
the L°° -norm of matrix valued functions is understood as taking the operator norm first and 
then the supremum over all (t, x) € C Then for any U £ C°°([0, T] x R")™ we have 

\\U\\Uc) < 2Te2^"^^)(|!C/|ii.(„„) + ||PC/||i.(£)). (8) 

a) Denote Qt := (0,i)xR" and 13 (t) := l + ||divyl(t, •) ~ B{t,-) - B*{t,-)\\^^.^^^y Then for any 
U G Ci([0,r],L2(R"))™nC°([0,T],iJi(R"))" the following estimate holds for allte [0,T]: 

t 

||C/(i,-)lli.(«„) < e^'^'^''(llC/(0,-)lli.(«™) + WPUWl.^n,))- (9) 



Proof. Applying the operator P to an arbitrary U G C^([0,T] x M"-)"^ we may write 

n 

{dtU, U) + J2{A^d^^ U, U) + {BU, U) = {PU, U) . (10) 

i=i 

A short calculation shows that 

n n 

2-Re{dtU, U) + 2Re^(^^a,^.C/, U) = dt\\U\\^ + ^ 9,^ {AW, U) - ((divv4)[/, U). 
Thus, taking two times the real part of (jlOl) we conclude 

n 

dt\\U\\^ + ^d^^{AiU,U) ~ {{'^^yA)U,U) + 2Re{BU,U) = 2Re{PU,U) . (11) 

To prove i), we rewrite ([TT|) as 

div(||^7||^ (AiC/, [/),..., (A„t/, U)) = ((divA)C/, C/) - {{B + B*)U, U) + 2Re(Pt/, t/) . 
Integrating over a partial lens £e ^ M"^^, the divergence theorem yields 

l{v%\\uf + Y^ vl^{AU, U)) dS = [{{divA -B- B*)U, U) dV + 2Re f {PU, U) dV . 



dCf 



i=i 



where vq is the unit normal vector field on dCs, assumed to point outwards with respect to Cq. 
Since vq = (— 1,0) , we conclude from inequality ([6]) that 



r " 1 

J i'^eWUf + ^ ^4(A^[/, U)) dS > -\\U{Q, •)||i.(«,) - ||C/(0, -n 



2 
dCe ^=^ 



Hence for all 6 G (0, 1] the term i||C/(6, ■)\\i^2,^^) is bounded by 

l|f^(0. ■)\\l-(no) + J ((div^ -B- B*)U, U) dV + 2Rc j {PU, U) dV . 

The terms on the right-hand side can be estimated with the help of the Cauchy-Schwarz inequality, 
leading to 

\\\U\\Une) ^ \\U{Or)\\Un,,) + ^imu\\l.^c.,) + \\PU\\l^^^^ 
and Lemmas 12.31 and 12.21 imply 

\mUc) < 2Te2T"(^)(|l[/(0, .)|li.(«„) + WPUWl^c,)- 
To prove ii), we integrate in pT|) over the spatial domain M", leading to 



i\\uit,-)\\l. 



J2 I dec, {AU, U) dVn + [{{divA -B- B*)U, U) dVn + 2Re j {PU, U) dVn 



After integration with respect to the time variable between and t we find 

t t 

||C/(t,-)lli2(R„) < ||f/(0,-)lli2(R„)+ /'/3(s)||f7(s,.)||i.(R„)ds+ f \\PU{s,-)\\l.^^„^ds (12) 



where we have used that Jg^„ dxj{A^{s,-)U{s,-),U{s,-))dV = for all j = l,...,n and for all 
s e [0,T] since a; i-^ U{s,x) belongs to iJ^(M")™ and the latter possesses C^(R")™ as a dense 
subspace. Employing Gronwall's lemma we turn (|12p into the desired estimate 



\\u{t, .)iii.(K„) < J'^'^'\\mo, .)iii.(K„) + iipc/|ii.(oo) 



n 



3 Generalized solutions to the Cauchy problem 

Having all necessary prerequisites at hand we draw our attention to the initial value problem ([TlH]) 
on the space-time domain fi^ := (0,T) x K". We will establish three statements of existence and 
uniqueness, each using different spaces of initial data and right-hand side. Working with a smaller 
space in this respect allows to relax the asymptotic conditions on the coefficient matrices A' and 
B. 

The formulation of the theorems requires some notions from Colombeau theory, we want to briefly 
review. A generalized function u S G{^) is called of L°°-type if it has a C°°-moderate rep- 
resentative (Me)e such that ||ue||^oom\ = 0{e~™) as e — > 0. It is called locally of logarithmic 
growth or locally log-type, if it has a C°°-moderate representative (Me)e such that for all K <Z<Z fi, 
||Me||^oo(^N = 0(log(l/e)) as e — )■ (cf. Definition 1.1 in [IS]). It is called of L°°-log-type if 
\\ue\\^^^^^ = 0(log(l/e)) as e ^ (cf. Definition 1.5.1 in [13]). A matrix A e Mm{G{fl)) is 
called Hermitian, if it has a Hermitian representative (^e)e, i.e. A^ is Hermitian for all e < Sq (cf. 
Lemma 4.3 in [53]). 

We call a partial differential operator P = dt + X]^=i A^dxj + B with Colombeau generalized 
coefficient matrices symmetric hyperbolic, if all matrices A^ are Hermitian and the entries of A^ 
and B are of L°°-iype together with all their derivatives, i.e. A^ ,B G Mm{GL°°(,^))- This ensures 
that there exists eo > and representatives {Al)^ and {B^}^ such that P^ = dt + X]?=i ^i'^xj + B^ 
is a classical symmetric hyperbolic operator for all e < Eq. The corresponding family of smooth 
solutions to the classical Cauchy problem for fixed e represents a candidate for the generalized 
solution. Yet some additional asymptotic growth conditions in e have to be imposed on the 
coeflacients to obtain a moderate family of solutions. In particular, certain log-type conditions on 
the coefficient matrices are essential in order to use a Gronwall-type argument in the proof (cf. 

The first theorem allows for the most general initial data and right-hand side, but requires the 
principal coefficients Al{t,x) to be bounded uniformly in e and {t,x) for large |a;|. 

Theorem 3.1. The initial value problem for a symmetric hyperbolic operator with Colombeau 
generalized coefficients, 

n 

dtU + '^A^dx^U + BU = F on ftr (13) 

U{0,x) = G{x), (14) 

has a unique solution U S G{Qt)"^, */ 

i) initial data G G 5(R")" and right-hand side F G g{[0,T] x K")™, 

ii) all spatial derivatives dx^A^ as well as the Hermitian part of B are locally of log-type, 
Hi) there exists Ra > such that \\Ai{t,x)\\ = 0(1) on (0,T) x {x e ]R"| |a::| > Ra} as e ^ 0. 



Proof. Wc pick Hermitian representatives (A^^ and representatives of B, F and G. There exists 
£o > such that Ve < Eqi the initial value problem 

n 
dtUe+Y.^edx,U, + B,U, = F, (15) 

J=l 

Ue\t=^ = Ge. (16) 

has a unique solution U^ G C°°{ViT)"^ (cf. Theorem 2.12 in 0). We claim that the equivalence 
class U = [{Ue)e] is the unique Colombeau solution. Hence we must show that the net {U^)e is 
moderate and that negligible variations of the coefficients and the data yield the same solution. 
Let K CC ftr and assume K C (0,r) x {x € M"| \x\ < Rk}- Choose i?i > max{RA, Rk) and 
^2 > ^1 + T{1 + 2y/nC). Then a standard lens C of thickness T, inner radius i?i and outer 
radius i?2 will contain K and satisfy inequality ([6]) by Lemma l2.4l Thus we may apply the energy 
estimate ([5]) e-wise and obtain 

||C/,||i.(^) < 2Te2r«^(^)(|lG,|li.(^„) + |lF,|li.(^)) = 0(e-™) (17) 

as £ — >■ for some tti G Nq, since the norms of the data grow only like some inverse power of s 
and as{C) ~ 1 + ||divyle ~ B^ — -^e IIl°°(£) ~ G(log(l/e)) as e — > by assumption ii). We attempt 
to show that all derivatives of U^ satisfy a similar estimate. For this purpose we introduce some 
convenient notations. For A e MmiC°°(nT)) and U G G°°(J7t)'" we define 



^\7^A 
T.^'A 



= d-mg{d,,A,...,d,^A) yic/ 

= diag(A,...,A)„. J:-U 



= (a,, [/,..., 9,„c/) 

= (c/,...,c/)^., 



T 



For example, VVl is a blockdiagonal matrix built from all spatial derivatives d"A of length \a\ = r. 
Similarly E'A is a blockdiagonal matrix whose blocks are just n^ copies of A itself. 

Claim 1. The vector VC/^ satisfies an equation of the form 

n 

dtWe + Y. Y^'^Ald^yu, + BIV^U, = Qr'S''C/e + VF, (18) 

i=i 

where B] = Yl^B, + (9,.Ai)i<,j<„, 4^+i = f.^Bl + (a,.EMJ)i<,,,<„ /or r > 1, and Ql~^ is 
a purely spatial partial differential operator of order r — 1 with coefficients depending linearly on 
spatial derivates of A^ and B^ up to order r. 



We present the case r = 1 in detail and proceed by induction. Differentiating (jl5|) with respect to 
Xk yields 

n n 

dtd-.,U, + ^ Ai9,^5,,C/e + Y, d^,A^,d^^U, + BA.Ue + d^,B,U, = d^,F, 

One would like to read this as an equation for dx^Us, but since the equations are coupled one has 
to consider the system 

n n 

There are no derivatives of U^ on the right-hand side and the only term that does not fit into our 
concept on the left-hand side can be rewritten in the following way, 

n 



so that V^Ue satisfies the system 

n 

where Bl = E^i?^ + (9a;iAi)i<ij<„ and Q° = — V^Sg. We proceed by induction with respect to 
the differentiation index r. Applying dx^ to the induction hypothesis (|T5| we find 

n n 

Just hke in the case r = 1 we try to write these k systems as one big system. It is easy to see 
that the right-hand side can be written as Q^^T.^^'^Ue + V+^F^ with Q^ a purely spatial partial 
differential operator of order r with coefficients depending linearly on spatial derivates of Ai and 
-Be up to order r + 1. Furthermore we can rewrite the lower-order terms on the left-hand side as 
B^+^V+if/e, which finally leads to 

n 

where indeed B^+^ = S^BJ + {dx,T,''A^)i<t,j<n- 

Since HS'yl^ll = 11^4^11^, the estimate ^ is also vahd for the symmetric hyperbolic operator 
P' = In-rndt + Ej=i S''-A^5a;^ +B''. By dll) in Lemma [23] i) we therefore have 

where 

<(£) := l + ||divSMe-i?,^-(i?nilL-(£) = 0(log(l/e)) as e ^ 0, 



since the Hermitian part of B^ can be constructed from the Hermitian part of B^ as well as first- 

2(£) 



order derivatives dx^Ai. From the fact that ||t^£||i2(£) — 0{e ™) and by iterative application of 



([TO]) for r = 1, 2, 3, ... we conclude that 

Vr e No 3m e No : || V''C/e||i2(£) = 0(e"") as e ^ 0. (20) 

It remains to show this asymptotic estimate for derivatives involving also the t-coordinate. 
Claim 2. All mixed derivatives d\^^Ue satisfy an equation of the form 

dlVUe = Rl^'^^'U, + dl-^VF, (21) 

where R^'^ is a linear partial differential operator of order r + I involving t-derivatives only up to 
order I — 1. Moreover, the coefficients of R^''' are linear combinations of spatial derivatives of A^ 
and B^ up to order r + 1 and time derivatives of A^ and Bg up to order I — 1. 

The case I — 1 follows immediately from ([T8)) by putting 

n 



since Ql ^ is a purely spatial operator of order r — 1 . Applying the operator dt to the induction 
hypothesis ((2T|) gives 

and i?J''+^I]''[/e :— dt{Rl''''S^Ue) is of course an operator of order r + I + 1 with time derivatives 
only up to order I. As an obvious implication of the Leibniz rule, its coefHcients are linear com- 
binations of derivatives of A^ and B^ with spatial derivatives of order r + 1 at most, since the 
coefficients of the operator Ql~^ depend (only) on spatial derivatives of ^^ and S^ up to order r+l. 

Successively making use of (|21l) for Z = 1, 2, 3, ... yields in combination with p^ that 

V;,r €No3meNo : ||(9tV''C/,|j^2(£) =0(e"™) as e -^ (22) 

and by the Sobolev embedding theorem on domains with locally Lipschitz boundary (cf. Theorem 
4.12 Part II in [T]) and the fact that K C C this implies 

VaeN^+i3meNo : ||9"[/e||i=e(;^) =0(e"™) as e ^ 0, (23) 

i.e. the class [(C/e)£] is moderate, since K CC ^t was arbitrary. For the uniqueness part we 
choose negligible nets (F^)^ and {Ge)e to represent right-hand side and initial data. From the 
energy estimates (flT)) . (TTOt . and equation ((2T|) it is then easy to see that the corresponding solution 
[(C/e)^] will also be negligible. By Theorem 1.2.3 in [13 it actually suffices to show the negligibility 
estimate for the zeroth derivative only, i.e. in terms of L^-estimates for all derivatives of C/^ with 
order < \{n + l)/2']. D 

If the 0(l)-condition on Al(t,x) for large |a;| is dropped, one cannot work with the same lens for 
all values of e anymore, but has to use an ^-indexed family of standard lenses. However, since 
the growth of the volumes of these lenses is under control, we can still keep many aspects of the 
solvability result by subjecting initial data and right-hand side to the stricter growth conditions 
of the space Gl'^^ ■ 

Theorem 3.2. Consider the alternative conditions 

i') initial data G G CJ^oc(M")™ and right-hand side F G Gl°^{^t)"^, 

ii') all spatial derivatives dx^A^ as well as the Hermitian part of B are of L°° -log-type. 

Then there exists a unique U G ^(fir)™ satisfying Equation Iil3\} and such that JJ^P holds in the 
following sense: U\t=o is equal to the image of G under the canonical map C/^oo(R")'" — >• t/(IR")™. 

Proof. The loss of condition iii) and the alternative version of i) have no effect on the applicability 
of Theorem 2.12 in [2 for fixed e < Eq. So we still get a solution candidate [(f/e)e] from the e-wise 
construction of a family (U^)^. After choosing a compact set K <Z<Z VLt, we again aim at building a 
standard lens around it such that ([6]) in Lemma [2^ is satisfied. Assuming that K C (0,r) x {x G 
R"! |x| < Rk}i we fix its thickness T and its inner radius Ri > Rk- Since we now have an 
£-dependent bound ||y4^||^oo < Ce~™, we choose a family of radii R2e ~ Ri+ T(l -1- 2-ynCe^™), 
thereby obtaining a family of standard lenses {C^)^ with outer radii R2e — 0(£^™) as e — > 0, all 
containing the compact set K. We put a^ = 1 + HdivA^ — B^ — ^eWL^^in )i i-^- the supremum 
taken over the whole domain fix. Employing estimate (|T7)) e-wise for the lens C^ leads to 



< 2re2^<(Vol„(Ho.)||Ge||i..(R„) + VoUi(£,)|li^,|l^^(^^)) = 0(£-™) 

for some ?ti G No, since both VolniHos) and Vol„+i(££) grow only like some inverse power of e as 
e —> and a'^ = 0(log(l/£)). Analogously using the higher order energy estimate ((19)) yields 

l|Vf/s|li.(^,)<2Te2^<'^(||V'^G£||i.(^„^)-M|Qr'S'^C^. + V'-+iF,||^.(^^)) 

< 4re2^<''(Vol„(Ho.)||V''G,||i^(R„) + ||Qr'S'-C/e||i.(^^) 

+ Vol„+i(/:,)||V''+i^,||^^(j,^)) = 0(e^™) as £ ^ 
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as the term HQ^ ^'E^Ue\\i^2(^ \ can be estimated via pulling out _L°°-norms of derivatives of the 

coefficients and a'/ := 1 + HdivEM^ - SJ' - {B^^)*\\^^ = 0(log(l/e)). With the help of ^ it is 
then easy to see that 

VZ,reNo3meNo : ||a^V''C/e||i2(£^) - 0(e"") (e -> 0). 

Since K C C^ for all e < Eq and by the Sobolev embedding theorem we conclude that for all 
a G Nq"*"^ there exists m e No such that |19"t/e||^oof^-) = 0(e~™) as £ — )■ 0. The uniqueness part 
is completely analogous to the corresponding part in the proof of Theorem 13.11 D 

Theorem 13.21 is applicable even for coefficients and data which are both associated to highly 
singular and periodic distributions. Denoting the delta distribution at y hy Sy — 6{- — y), it 
is possible to consider, e.g., principal coefficients A^{t,x) « X^ksn™ '^''«(^)^''(*) ^^'^ initial data 
G{x) ~ X^KgN™ '^9«(^) '^ith real numbers /,g > 0, representing m-dimensional lattices of Dirac 
measures with lattice constants l/l and l/q, respectively. 



Remark 3.3. The conditions in Theorem \8.S\ allow for infinite propagation speed near spatial 
infinity as e -^ Q. In general this may cause non-uniqueness of solutions, see Example 17.1 in 
\27l. Using the space Ql"^ for initial data and right-hand side avoids nonuniqueness, yet null 
"solutions" with non-vanishing initial data still exist. In fact, any initial data G in the kernel of 
the canonical map ^l=(R")"' -^ ^(M")™ yield [/ = oQ 

Initial data and right-hand side decaying at spatial infinity (|a;| — >■ oo) make it possible to relax 
the conditions on A' and B even a bit further. More precisely, the required asymptotic behavior 
of Ai and i?e + B* can be made less restrictive with respect to the time variable. We use mixed 
norms 11^11^1. =0(^2^) '■= Jq \\A{s, •)I1l=°(R")'^* fo^ ^^Y ^ ^ Mjn{C^ (Ht)) ■ We say that an element 
A e Mm{GL°^ (^t)) is of L-'^'°°-log-type, if it has a representative {A^)^ such that H^ell^i.^offj ) = 
0(log(l/e)) as e ^- 0. A similar norm was introduced in Definition 2.1 in [S]. 

Theorem 3.4. In the initial value problem IIS^H^ assume that 

i") initial data G £ 5l2(R")" and right-hand side F G 5i2(f]T)'", 

ii") all dxiA^ as well as the Hermitian part of B are of L^'°° -log-type. 

Then there exists a unique solution U G ^^2(1}^)™ to the initial value problem hlS^H^ . 

Proof. Fixing Hermitian representatives of A^ and representatives oi B, F and G, we may use 
Theorem 2.6 in [5] to provide solutions U^ G C°°([0,T],i7°°(R"))" to the classical initial value 
problem for each e < £o- To show moderateness, we plug U^ into the energy estimate (|9]) and find 



|C^e|li2(n,)<r sup ||;7e(t,-)lli2(K„) 
0<t<T 



<Tef'^{\\Ge\\U^., + m\\UnJ=0{e-n, (24) 



where Pe := ||divAe - B^ - i3*||^i,oo(n.j,) = 0(log(l/e)). In Claim 1 in the proof of Theorem [XTI 
it has been shown that the vector V^Us satisfies an equation of the form 

n 

dtV^U, + Y, J^^A^.d^yU, + B^VC/e = Qr^S'^t/, + V^e (25) 



^Consider the scalar equation dtu H dxU = with u\t—o = i^, where 95 G I5(R) and ¥3(0) = 1. Then the 

unique Colombeau solution is the class {{^(x i))e] = in C7((0, T) X R) since for all K CC (0, T) X R we have 

supp(i/p(a; — ii) n X = for e small enough. 
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where Ql ^ is a purely spatial partial differential operator of order r — 1, SM.i are simply block- 
diagonal matrices consisting of A^-blocks and B^ depends solely on B^ and spatial derivatives 
d^.Ai. Thus, plugging V'T/e into dH) yields 

^sup^ ||V'^C/,(t, .)|li.(M„) < e^= (||V'^G,||i.(K„) + IIQr'S'^C^. + V'^F,||^.(^^)), (26) 

where ^l = l + ||divE''^£ - B^ - (-B^)*||^i.oo(fi^). Iteratively applying ([26]) shows that for all r e N 
there exists m e No such that ||V''C/e||^2(j2^) ~ 0(e^™). Finally, successively employing equation 
(PT|) for I — 1,2,3,... one finds that the asymptotic growth of \\dt'^^Ue\\j^2rQ\, \\df\7^Ue\\i^2rQ\, 
\\dfV^Us\\j^2/^2 )' ■■■ ^^ ^^^'~' moderate. Altogether we therefore have 

Va e N[J+i 3m e No : ||a"C/,||i.(j,^) = 0(£-™) as e ^ 



and hence [(C/e)e] € Gl^{^t)- To show uniqueness, we assume negligible data (F^)^ e A/'l2(^t) 
and {Ge)e G A/l2(]R"). The energy estimates used to prove moderateness then immediately imply 
negligibility of the corresponding solution [{Ue)e]- D 

Note that in Theorem 13.41 thanks to the L^'°°-norms in condition ii") no logarithmic scaling of 
the moUifier is required to model a lower order coefficient of the form B(t, x) ~ 6(t)B{x) where S 
represents the delta distribution B is bounded. 

Remark 3.5. For the existence and uniqueness results presented in this section, there exist ver- 
sions which are also global in time. In correspondence with Theorem \S.1[ given coefficients A' 
and B m M™(5i,oo(R"+i)), data F e 5(M"+i)'" and G e ^(M")'", one obtains a global solution 
U € tJ(K"+i)"' if all dxiA^ as well as the Hermitian part of B are locally log-type and A^ = 0{1) 
as e — >■ outside a cylinder R x -B_r^ . Analogously extending the respective asymptotic growth 
conditions from Ht toM"^^ yields "global in time" -variants of Theorems W^ and \3^\ Concerning 
"global in time" -solutions, the requirements on the coefficients with regard to the dependence on t 
can be somewhat relaxed. It suffices to demand the estimates in the assumptions of the theorems 
only locally in time, i.e. for all I CC M. A similar observation was made in Remark 1.5.3 in flSf . 

4 Regularity of the generalized solutions and distributional 
limits 

To justify the term "generalized solution" , compatibility with classical smooth and distributional 
solutions should be investigated. When the coefficients are smooth, the unique Colombeau solution 
should be equal to the respective classical solution in a certain sense. As in 16.^ and |23J, the fol- 
lowing proposition establishes compatibility with the classical results for smooth and distributional 
data. 

Proposition 4.1. (i) In Theorem \3.1[ additionally assume that A^ and B have components in 
C^i^T)- If F e C°°(f7T)" and G e C°°(IR")" then the generalized solution U G ^(f^r)™ is 
equal to the classical smooth solution. 

(li) Suppose thatAi and B m Theorem\^are smooth. For s (^M. let Fq e L2([0,r],7J"(M"))™ and 
Go e H''{W)"' and denote by Uq the unique distributional solution to (EP-|0i in C([0, T], iJ''(R"))'". 
Define generalized data by F := [iFs)e] G Sl^COt)" and G := [(Ge)^] G ^l2(R")'", where F^ and 
Ge are moderate regularizations such that F;, — >■ Fq in Z/^([0,r],.ff*(]R"))'" and G^ — > Go in 
if''(R")'" as £ — >• 0. If U = [{Ue)e\ is the corresponding generalized solution in gi^2{QT)"\ then 
Ue -^ Uo in G([0,r],iJ''(R"))". 

Proof, (i) Since we may choose the constant nets {F)^ and (G)e as representatives of the classes 
of F and G in G, we obtain the classical smooth solution to problem ([TJ2]) as a representative of 
the unique Colombeau solution. 
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(ii) By Theorem 2.6 in 3] we have for all t (E [0,T] 

||f/e(i,.)-f/o(t,-)llL(E") 

t 

< C{\\G, - Go||?,.(j,„) + I \\F,{s, •) - Fo(s, ■)\\l.^^,.^ds), 



where the constant C does not depend on e, since all coefficient matrices A^ and B are assumed 
have components in C°°(17t). Taking the supremum over all t € [0, T] and letting e —>■ gives the 
convergence in C([0, T], i7^(M"))™. D 

We may interpret Proposition 14. 1 f ii) as a statement on the regularity of the generalized solution, 
measured in terms of iy-norms of the associated distribution. Yet this concept of regularity is 
restricted to situations where distributional limits exist and therefore not applicable if initial data 
or right-hand side are not associated to any distribution. Intrinsic regularity theory in Colombeau 
algebras is based on the subalgebra Q°°{n) of regular generalized functions in G{^) and has been 
investigated in the context of hyperbolic partial differential equations in [31] [ini [2H1 UHl IH] • In the 
study of intrinsic regularity of generalized solutions to partial differential equations, the notion 
of slow scale nets was introduced in [20 and has proven to be essential in many circumstances 
(cf. [22 [ini E])- A net {rs)s of complex numbers is said to be of slow scale if \rs\P — 0{l/e) as 
e — ^ for all p > 0. As in J16j we call a net {se)e of complex numbers a slow-scale log-type net 
if there is a slow scale net (r^)^ of real numbers, r^ > 1, such that \s^\ — 0{\og{r^)) as e — > 0. 
Generalized functions satisfying the moderateness estimates with slow scale nets in place of the 
inverse powers e~^ are called slow scale regular. They represent a different notion of regularity 
than regular generalized functions. 

Proposition 4.2. In Theorem \3.1[ assume all coefficients A' and B to be slow scale regular 
generalized functions. In addition suppose that all log-type conditions are replaced by slow-scale 
log-type estimates. Then F € g°°{[0,T] x R")™ and G G 5°°(M")™ implies U G ^""(^Jt). 

Proof. Fix a standard lens C of thickness T with initial surface "Ho- Then there exists M E No 
such that for all r, Z, we have 

||alVF,||^.(^) = 0{e-^') as well as || VGe||^.(^^) = 0(£-*^). 

Since all coefficient depending factors in the L^-estimates (|T7)) and (IT^ are of slow scale, for each 
a e Nq we obtain a certain slow-scale net (re)^ of positive real numbers such that \\d^Ue\\i^2(£\ = 
0(rge^*^). Finally, by (|2T|) we only pick up slow-scale factors with each time derivative applied 
to d"U^, and so we have 

for all I E No and a e Nq . This proves the assertion. D 

The intrinsic regularity property holds also for the generalized solutions obtained from Theorems 
13.21 and 13.41 respectively, if all log-type conditions are replaced by slow-scale log-type estimates 
and all coefficients are slow scale regular generalized functions. Note that all these conditions are 
automatically satisfied when the coefficients are smooth (like in Proposition 14. II (i)). 
Despite its t/°° -regularity, the solution in Proposition l4 . 2 I mav not be associated to any distribution. 
In the remaining part of this section we want to investigate distributional limits of the generalized 
solution when the coefficients are non-smooth. To this end we consider the generalized Cauchy 
problem (J13m4p with data and coefficient matrices as regularizations. For convenience of the reader 
we first collect some basic properties of regularizations of W^'^''P-functions obtained by convolution 
with a moUifier. 

Lemma 4.3. Let p e I?(K"), J p{x)dx = 1, supp(p) C i?i(0), and put Pe{x) — e^'^p{x/e). Given 
u G IF'^'P(R"), k £ No, 1 < p < 00 we put Ue :— u* p^. Then we have 

13 



i) [("s)e] e eL~(ffi") and |19"we|lty.,p(R„) = 0(£-l"l), 

iij ||we — u\\y^,k,q,^r^\ — i' as e -> /or u e P^'^''^(]R") mi/i 1 < (7 < oo, and 

Hi) if u is bounded and uniformly continuous, then \\u^ — u\\i^somn) — > as e — >• 0. In particular, 
for u e W^'''°°(K") we have \\ue - u||vi/fc-i,oo(R,i) ^ as e ^ 0. 

Proof. Here, all norms are taken on the domain R". (i) If w e VF*'''^(R") and a G Nq we may esti- 
mate ||9"Ue||^oo = IK^^pe) * u||^oo < ||5"p£||^p'||w||^p with l/p+ 1/p' = 1 by Young's inequality 
and we have d"pe = e-l"l(5"p)e which yields Wd^u.W^^ = 0(e-l"l), hence [(u,),] £ ^loo(R"). 
Moreover we find 

lia"«,ll^.,p = max||a"+VlL. < max ||(a>,) *a^u||^, < "^°^"^Jj|""^'•^^ 

For (ii) and (iii) we refer to [Sj Theorem 8.14]. D 

The following proposition is concerned with the generalized Cauchy problem (fT5|) - dTil) when both 
data and coefficient matrices are obtained from convolution regularizations of non-smooth coeffi- 
cient matrices Aq and Bq, right-hand side Fq and initial data Gq. 

Proposition 4.4. In the initial value problem (OST^Pi assume that the entries of Aq and Bq are 
bounded and Lipschitz continuous with respect to X, Gq e i/^(R")™, and Fq £ L^([0,r], iJ^(R"))'". 
Define corresponding generalized coefficients and data satisfying the assumptions in Theorem \3.4\ 
H Then there exists Uq £ C([0, T], ^^(K"))™ n H^iflT)"" such that U^ -^ Uq in the norm 
L°°([0, T], i^(R")) for every representative {Ue)e of the unique Colombeau solution U £ Qi2(flT)™'- 

Proof. First observe that ||9jC4||^2m ) — 0{1) as e — > for all j. For the spatial derivatives 
dx U^ this follows from the estimates (IM)) and (0^ by virtue of the uniform boundedness of 
/3l = 1 + ||div "E^Ag — Bl — {Bl)*\\^x^^,^ -, as well as {Ge)e and (Fg)^ in the respective norms (see 
Lemma [4. 3p . The L^-norm of the time derivative dtUe can then be estimated directly by means 
of the equation. We show that the e-wisc constructed family of solutions to the Cauchy problem 
(HaHl), {Ue)e, is a Cauchy net in C{[Q,T],L'^{W))'^. For indices < e < e smaU enough we 
obtain from 

sup ||f/,(t, •) - U,{t, •)IL.(K„) < Ct{\\G, - G,|L.(K„) + ||F, - P,C/,|L.(o^)), 

Q<t<T 

where Ct is independent of e and Pe = dt + Ej=i ^i^x^ + B, = dt + ^"=1 {M ~ ^ + ^D^x, + 
B^ — Bg + Bg. Thus, the energy inequality yields 

\\Ue - C^£|Il=-([0,T],L2(R")) < CryllGe - Ge||^2(R„) + \\F^ - i^e||2,2(j2r) 

n 

where we omitted the domains in the norms on the right-hand side (all norms are taken over the 
strip Qt or R" respectively). The convergence properties of the nets {Ge)e, {Fe)£, i^De, {Be)e 
and the uniform boundedness of llc^^jC^ell r2/Q -) imply that {Ue)e is a Cauchy net in the norm 

L°°([0, T], L2(R")), thus there exists Uq £ C([0, T], L2(R"))™ such that \\U, - t/o|lL~([o,T],L2(K")) ^ 
as £ —^ 0. Since \\djUe\\i^2u-i \ < C, we have for any ip £ I?(J7t)™, 

Hence djUo £ L^iW')"", thus Uq £ C{[0,T], L^{W'))"' n iJi(^T)". □ 



^The conditions in Thcoreni l3.4l are automatically satisfied when regularizing via convolution with a nioUifier as 
in Lemma 14. 31 
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The conditions on the coefficient matrices in Proposition 03] are typical assumptions in results on 
weak solutions. For a comparison of solution concepts for linear first-order hyperbolic differential 
equations with non-smooth coefficients we refer to [14 . 

Corollary 4.5. The initial value problem (MW with Aq, Bq, Gq and Fq as in Proposition \4-4\ 
has a unique weak solution Uq G C([0, T], L^(R"))™ n H^{VIt)"^, equal to the distributional limit 
of the generalized solution. 

Proof. We show that the distributional limit of the generalized solution U = [(C/e)^] is the unique 
weak solution. First observe that the limit Uq of U^ is continuous in time, so it satisfies the initial 
condition, i.e. C/o|t=o — Go- Moreover we have C/g -^ Uq in L'^{nT) and thus also Bi^U,. -^ BqUq 
in L^(J7t)™. The proof of Proposition 14.41 shows that |jC4||//i(si ) is uniformly bounded. By 
the weak compactness theorem there exists a weakly* convergent subsequence (t/j_)/cgN with 

limit Uo e H^inr)"" (cf. Theorem 6.64 in ^). Thus dtU_L. -^ dfUo weakly* in L^iVtT)"" and 

Uo = Uoe C([0,T],i2(M"))™ni/i(f7T)™. 

We have \{d,,U^,iA'^ - Al)4,)\ < \\d^^U^\\ ||(A^_^ - Al)^P\\ for any tA e L^^rr, 

hence {A'^d^^U^,^) = {d^^^M'^ - K)i') + (5., t/^ , ^'oV') converges to (9,^.[/o, A^oV-)- 
Therefore we obtain 

k 



lim 



1 (^{dtU„ ^) + Y,{Aid,,U„ ^) + {B,Ue, ^)) = lim(F„ .A) = (Fo, tA). 



To prove uniqueness, suppose Vq € C([0, T],L'^{W^))'^C^H^{VtTT is a solution such that PVq = Fo 
and Fo|f=o = Gq. We may regularize this solution so that Ve -> Vo in C([0, T],i^(R"))'" and in 
H^inr)"^. Then Ve\t=o ^ Go in L2(R«)m and necessarily lim (^tK + Ei=i M^x^V^ + B^V.) = 

Fo in ^^($7^)™ irrespective of the regularizations chosen. Applying the basic energy estimate ([9]) 
to the difference Ue — V^ yields 

sup \\Ue{t, •) - V,{t, OIL.fK,.) < Ct{\\G, - Fe||^.(R„) + ||F, - FV;||^.(,,^)) ^ 

0<t<T 

as £ -> and therefore Uq = Vq in C([0, T], F2(m«))™. q 

Note that the statement in Corollary 14.51 is not just a compatibility result, in fact we directly 
obtain a unique weak solution to the initial value problem with non-smooth coefficients merely 
by studying properties of the generalized solution when the coefficients and data are regularized 
distributions. In case of smooth coefficients, the methods applied in the proofs of Proposition 14.41 
and Corollary 14.51 would lead to the corresponding classical existence and uniqueness results as 
well (see Proposition 14. ip . 

Remark 4.6. (a) In order to get more information on the regularity of the distributional shadow 
Uq and its dependence on the regularity of the coefficients Aq, Bq and data Fq, Gq, it is essential 
to have precise estimates on the speed of convergence of the regularized objects (cf. the notion 
of "strong association" in JSTI). For example, rapid convergence in the principal part guarantees 
the existence of a distributional shadow under regularity assumptions on the initial data Gq and 
right-hand side Fq which are weaker than those in Proposition \4.4\ To be more precise, we have the 
following statement: In Theorem \8.4\ let all Aq be Lipschitz continuous and let Bq be uniformly 
continuous and bounded. Given Gq G F^(R")™ and Fq g L'^{flT)™, define generalized coefficients 
A, B and data F, G such that as s — > 
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Then the corresponding generalized solution U = [{Ue)e] still has a distributional shadow Uq G 

L'^{Q,tY"' such that Ue — >■ Uq in L^^flx)"^ for any of its representatives (Ue)e- 

The proof strategy is the same as in Proposition \4-4\ showing that (f/e)e is a Cauchy net in 

L'^^Ht), the main difference now being that dx U^ is not bounded, but kept under control by the 

factors stemming from the coefficients. 

(b) To illustrate that the conditions discussed in (a) are realistic, we gather the following estimates 

on the convergence speed of convolution regularizations: Let p G iS(M") with J p{x)dx — 1 and 

J x°'p{x)dx = for all |a| < ?7i G N, s G R and put p^ := (T^"p(^) where CTg — > as £ — > 0. Then 

we have 

i) ue VF™'°°(R") => \\pe * w - w|Il-(R") = 0(cr™) as e^O, 

a) u G i^(M") =^ llpe * u — w||^a(R„) = 0(cr™) as £ — > for all s < —^ - 1 - m. 

Hi) u G iJ^(M") =^ ||pe *u — u||^.wjj„> — ^ and ||pe * u — wH^s-™™,.-) = 0((t™) as £ ^- 0. 

Here, i) is shown by Taylor expansion of u, ii) can be proved by considering the action on a test 
function and Taylor expansion of the latter, and Hi) follows from the definition of the Sobolev 
norms in terms of the Fourier transform. 
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